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ABSTRACT 



^ I We present model calculations of ionization fractions for elements in the 

OO I warm (T ~ 10*^ K), low-density photoionized interstellar medium (WIM) of the 

^ ' Milky Way. We model the WIM as a combination of overlapping low-excitation 

^ ■ H II regions having n{H^)/n{H) ^ 0.8. Our adopted standard model 



incorporates an intrinsic elemental abundance pattern similar to that found for 
warm neutral clouds in the Galaxy and includes the effects of interstellar dust 



S ■ grains. The radiation field is characterized by an ionizing spectrum of a star 



O ■ 

^ ■ with Teff ^ 35, 000 K and an ionization parameter log(g) ~ —4.0. The emergent 

emission line strengths are in agreement with the observed ratios of [S II] /Ho;, 
^ - [N II] /Ha, [S II] /[N II], [O I] /Ha, [O III] /Ha, and He I/Ha in the Galactic 

' WIM. Although the forbidden emission-line intensities depend strongly upon 

the input model parameters, the ionization fractions of the 20 elements studied 
in this work are robust over a wide range of physical conditions considered in 
the models. These ionization fractions have direct relevance to absorption-line 
determinations of the elemental abundances in the warm neutral and ionized 
gases in the Milky Way and other late-type galaxies. We demonstrate a method 
for estimating the WIM contributions to the observed column densities of 
singly and doubly ionized atoms used to derive abundances in the warm neutral 
gas. We apply this approach to study the gas-phase abundances of the warm 
interstellar clouds toward the halo star HD 93521. 



Subject headings: ISM: abundances - ISM: atoms - ISM: H II regions - Galaxy: 
abundances - Galaxy: general - radiative transfer 
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1. Introduction 

In elemental abundance studies of the warm, neutral interstellar medium (WNM), it is 
generally assumed that ionization corrections for photoionized gas are negligible, or at least 
modest. In practice, this assumption usually takes the form N{X)/N{H) ^ N[X'^)/N{H^), 
where X' is the dominant ionization stage of element "X" in the neutral gas, even though 
X* may also be found in ionized regions (see Savage & Sembach 1996). Ignoring ionized gas 
contributions to N(X'^) is often justified when studying cold atomic and molecular clouds, 
but this simplification introduces additional uncertainties into studies of the interstellar 
medium (ISM) in the immediate vicinity of early-type stars or along extended, low-density 
sight lines. In these latter cases, ionized material can represent a significant fraction of 
the total amount of interstellar gas observed and can contribute to the column densities of 
singly and doubly ionized atoms. 

The warm (T~10^ K), ionized interstellar medium (WIM) outside classical H II regions 
is a fundamental gas-phase constituent of the Milky Way and other late-type spiral galaxies. 
In the Milky Way, the WIM is detected through faint optical emission-line radiation and 
pulsar dispersion measures, which together reveal that the ionized gas has an exponential 
scale height of ~ 1 kpc, a volume fiUing factor of ^ 20%, and an average electron density 
of ~ 0.1 cm~^ (Reynolds 1993). Though not as easily observed as the H I layer revealed 
through 21cm emission, the WIM has a mean mass surface density (~ 1.6 M© pc~^) 
nearly one-third that of the H I distribution (Reynolds 1993; Dickey & Lockman 1990). 
Pervasive warm ionized gas distributions have been observed in other spiral galaxies besides 
the Milky Way to varying extents (Dettmar 1990; Rand, Kulkarni, & Hester 1990, 1992; 
Rand 1996; Hoopes, Walterbos, & Rand 1999). In all cases it seems clear that while there 
may be a variety of ionization sources and mechanisms contributing to the WIM emission, 
photoionization by OB-stars is the only source of energy sufficient to maintain the bulk 
ionization properties of the WIM (which in the Milky Way requires ^ 10~^ erg cm~^ s"-*^, 
or ~ 15% of the ionizing radiation output from OB-stars). 

There are two recent observational motivations for modeling the WIM and its impact 
on elemental abundance studies of the ISM. First, there is an increasing database of 
emission line measurements of the WIM. In particular, the Wisconsin Ha Mapper (WHAM, 
Reynolds et al. 1998b) is providing wide sky coverage of Ha and other optical emission hnes 
in the diffuse gas. Second, the Far Ultraviolet Spectroscopic Explorer (FUSE) will soon 
produce high spectral resolution absorption-line measurements in the 905-1195 A bandpass, 
which contains many diagnostics of photoionized gases. Together with absorption-line data 
from the Hubble Space Telescope Imaging Spectrograph (HST/STIS), FUSE will enable 
comprehensive studies of the elemental abundances in the WNM and WIM. 
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In this paper, we estimate ionization fractions for the WIM based upon existing 
knowledge of its properties. These ionization fractions can be used to derive appropriate 
corrections for ionized gas contributions to the measured amounts of singly and doubly 
ionized species found in both the WNM and WIM. We begin in §2 by summarizing the 
current observational constraints on the WIM. In §3 we describe our model of the WIM, 
its boundary conditions, the parameter space explored, and the sensitivity of the predicted 
emission line ratios and ionization fractions to the input parameters. In §4 we provide 
a short description of the applicability of the model results to absorption-line studies of 
elemental abundances in the ISM. Section 5 contains a brief discussion of these results and 
a few concluding remarks. Ionization fractions for radially-averaged sight lines through low 
density H II regions can be found in Appendix A. 

2. Emission-Line Properties of the Warm Ionized Medium 

The primary optical emission hues used to characterize the WIM are Hq;A6563, 
[S II] A6716, [N II] A6583, [O I] A6300, [O III] A5007, and He I A5876. We compile a set 
of measurements of these species appearing in the literature in Table 1, where we list 
the line intensities normalized to the observed Ha intensity. Compared to classical H II 
regions, the WIM has high [N II] /Ha and [S II] /Ha and low [O III]/H« ratios (Reynolds 
1985b; Osterbrock 1989). A subset of the listed ratios has been modeled and discussed by 
Domgorgen & Mathis (1994), who found that dilute stellar radiation leaking through the 
ISM provides an adequate source of ionization for the WIM. This paradigm is supported 
by models of the radiative transfer of ionizing photons through the Galaxy subject to the 
presence of interstellar clouds (Miller & Cox 1993; Dove & ShuU 1994; Bland- Hawthorn 
& Maloney 1999). Additional secondary ionizing sources may be present; in at least one 
case (Ogden & Reynolds 1985), the emission from a filament within the diffuse background 
emission appears to be due to coUisional ionization by a weak shock. However, the large 
[N II] /Ha ratio is atypical of the more general value of 0.4-0.6 for the emission observed in 
other directions (see Table 1). 

The emission-line ratios for gas in the Local and Perseus spiral arms change significantly 
as a function of the observed Ha intensity (Haffner, Reynolds, & Tufte 1999). Ratios of 
[S Il]/Ha and [N Il]/Ha in excess of unity are seen in some directions where 1^^ < 0.5 R. 
The ranges of [S II] /Ha and [N 11] /Ha listed in Table 1 for the Perseus arm encompass 
a majority of the observed values and are appropriate for heights z ^ 1.0 kpc above the 
Perseus arm. Haffner et al. (1999) argue that the [N 11]/Hq; ratio probes the temperature 
of the gas, while the [S II] /[N II] ratio traces the ionization state of sulfur. (The [S II] /Ha 
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ratio is therefore a probe of both physical properties.) They interpret the large, systematic 
increase in [N II] /Ha with height above the plane as an indicator of increasing temperature, 
with an as yet unidentified source of heating at large distances from the Galactic midplane. 
Extragalactic background radiation is an unlikely source, since the ionizing flux out to 
about 100 kpc from the Galactic disk is dominated by light from the Galaxy (Maloney 
& Bland-Hawthorn 1999). Interestingly, even though the [S II]/Hq; and [N II] /Ha ratios 
change by more than a factor of three with height above the plane in the Perseus arm, the 
ratio [S II] /[N II] ?a 0.5 — 0.65 is relatively constant and shows virtually no dependence on 
height above the Galactic midplane. Thus, in the interpretation presented by Haffner et al., 
the increase in temperature with height above the plane is not matched by corresponding 
changes in the ionization state of the gas. 

To date, the strongest constraints on the ionization of the WIM come from [O I] A6300, 
[O III] A5007, and He I A5876 emission-line observations (see Domgorgen & Mathis 1994). 
A complete set of optical emission-line measurements is not yet available even for a single 
direction; we make the simplifying assumption that the range in values and upper limits 
listed in Table 1 are representative of the WIM as a whole. However, we also recognize 
that a range of physical conditions is possible and allow the model parameters to vary 
accordingly. Additional emission-line data in the near future will hopefully lead to a better 
description of the intensities of the less well-studied species. 

3. Model Description 

To study the physical conditions of the Galactic WIM, we use the ionization equilibrium 
code CLOUDY (v90.04; described by Ferland et al. 1998 and Ferland 1996) to derive 
the temperature and ionization structure of spherically-symmetric nebulae surrounding 
single ionizing sources. We model the WIM as a combination of low-density, low-excitation 
H II regions. Our goal is to estimate the ionization fractions, x{X^) = N{X^)/N{X), 
of several ionization stages, i, of the most important elements, X, for studying the 
gas-phase abundances of the Galactic ISM. In doing so we endeavor to match the observed 
emission-line observations of the WIM as closely as possible. 

3.1. Approach and Assumptions 

Following Domgorgen & Mathis (1994), we model the WIM of the Galaxy as a 
combination of low-density H II regions. Our models calculate the temperature and 
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ionization structure of the model H II region from a distance of 0.3 pc from the exciting 
source to the point where the fraction of neutral hydrogen, x{H°), at the edge of the 
nebula becomes larger than a critical value, x^dge- Domgorgen & Mathis (1994) used models 
with Xedge = 0.10 to approximate fully-ionized material and those with Xgdge = 0.95 to 
simulate ionized interfaces of neutral clouds or other regions with a significant fraction of 
neutral hydrogen. For continuity with their work, we consider the same x^dge values in 
this study. In some situations, x^dge values as high as 1.00 may be encountered as the 
sight hne passes into a fully neutral cloud, but the model results for Xedge — 0.95 closely 
approximate those for Xgdge — 1-00. The emission line intensities of the gas in these two 
cases are indistinguishable, and the ionization fractions vary by less than ~ 10% for most 
ions. In this study we also consider a hybrid gas mixture described by a linear combination 
of equal amounts of Xgdge — 0.10 gas and Xedge — 0.95 gas. 

The characteristics of a low-density photoionized nebula are uniquely determined by 
three input parameters, assuming a given set of atomic data. These parameters are the 
set of gas-phase abundances adopted for elements heavier than hydrogen, the shape of 
the ionizing spectrum, and the ionization parameter, q, which is a measure of the ratio of 
ionizing photons to particles in the nebula (see Domgorgen & Mathis 1994 and Howk & 
Savage 1999). The inclusion of dust grain heating and cooling can also affect the model 
results, but these effects are small in the models considered (see §3.3). 

We consider three sets of abundances in our models: solar system, B-star, and warm 
neutral interstellar medium values, all of which are summarized in Table 2. The abundance 
of He is assumed to be one-tenth the abundance of H in all three cases. The solar system 
abundances are meteoritic values from Anders & Grevesse (1989) except for C, N, and 
O, which are photospheric values from Grevesse & Noels (1993). The B-star abundances, 
which may represent a better "cosmic" reference system than solar system abundances (see, 
e.g.. Savage & Sembach 1996 and Mathis 1996), are averages of the Kihan-Montenbruck, 
Gehren, & Nisscn (1994) determinations. We assume that the B-star abundances of 
elements lacking observed values are equal to the solar values scaled downwards by 0.2 
dex. We also consider a WNM-like abundance pattern to account for the possible effects of 
elemental incorporation into dust grains in the ionized interstellar medium. The elemental 
depletion pattern inferred for the WIM is comparable to that observed in the WNM (Howk 
& Savage 1999; see also Lagache et al. 1999). We adopt WNM gas-phase abundances 
equal to those observed in the well-studied warm diffuse clouds toward the low-halo star 
fjL Columbae (Howk, Savage, & Fabian 1999). The observed abundances in these clouds 
exhibit the imprint of elemental incorporation into dust grains. We adopt C, N, and O 
abundances for the WNM from Sofia et al. (1997), Meyer et al. (1997), and Meyer et al. 
(1998), respectively. 
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Our models are unable to match the [S n]/[N II] ratios observed by Haffner et al. 
(1999) if we adopt a solar abundance for sulfur. The predicted vahics of [S II]/[N II] are too 
high by roughly a factor of 1.5. We have explored possible explanations for high [S Il]/[N II] 
ratios within the context of our model, but find that the most reasonable explanation is a 
sub-solar reference abundance. (Haffner et al. (1999) reached a similar conclusion.) The 
/i Columbae clouds exhibit nearly solar gas-phase abundances of the non-depleted elements 
S, Zn, and P, but this sight line may be unusual in this respect (Howk et al. 1999). 
Therefore, for the elements S, Zn, P, Ar, and Na we adopt WNM reference abundances 
equal to the solar system values scaled downwards by 0.2 dex, in agreement with our current 
understanding of the abundances in the local ISM (Mathis 1996; Snow & Witt 1996; Meyer 
et al. 1998). This is essentially equivalent to adopting a B-star reference abundance with 
minimal incorporation of these elements into interstellar dust grains. 

All of the H II region models calculated with the adopted WNM abundances include a 
standard mixture of interstellar graphite and silicate grains (see Ferland 1996 and Baldwin 
et al. 1991). Dust grains in the model nebulae affect the thermal structure of the medium 
through photoelectric heating. The grain opacity also alters the shape of the ionizing 
spectrum slightly. The incorporation of potential coolants into grains is an essential 
component of the model and is required to match the observed emission-line strengths. 

We use ATLAS line-blanketed LTE stellar atmosphere models (Kurucz 1991) as input 
spectra to the CLOUDY ionization code. The input spectral shape is varied by changing 
the effective temperature of the stellar atmosphere of the central ionizing source. We 
consider model H II regions ionized by central stars with stellar effective temperatures in 
the range 31, 000 < T^ff < 41, 000 K. 

The third input to our models, the ionization parameter, is a dimensionless measure 
of the relative particle and ionizing photon densities. Throughout this work we write the 
ionization parameter, q, in the form 



where L50 is the stellar ionizing luminosity in units of 10^*^ photons s~^ and / is the 
volume fiUing factor of the ionized gas (Domgorgen & Mathis 1994). This definition of 
the ionization parameter, q, is related to the "volume averaged" ionization parameter. 



q = nnfL^o, 



(1) 



U — |(n-y/ne), through 
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where ob is the recombination coefficient of H to levels n>2, and is the number density 
of ionizing photons. 

Bright high-density H II regions typically have values of log(g) ^ —1.0. We calculate 
models over the range —4.0 ^ log(g) <^ —3.0, which is more appropriate for the low-density 
WIM (Domgorgen & Mathis 1994). For a nebula surrounding a normal 09 V star with 
Teff = 35,000 K, total luminosity logL^,/LQ = 5.0, and ionizing photon flux logLso ~ —1.6 
(Vacca, Garmany, & Shull 1996), a value log(g) = —4.0 corresponds to n^P ~ 4 x 10^'^ 
cm~^. Estimates for the density of the WIM typically yield fue ~ 0.08 and / ^ 0.20 
(Reynolds 1993). Note that the ionization parameter rather than the density determines 
the properties of the model nebulae; for a given value of q in Eq. (2), higher values of 
can be achieved provided there is a proportional change in n^. 

We extract from our models the volume-averaged intensities of various emission 
lines relative to Ha, as well as the ionization fractions, x(X*), of elements important for 
absorption-line studies of the low-density Galactic ISMQ. We use volume-weighted averages 
rather than radial averages since a random, extended sight line through the WIM is 
likely to sample gas found within a variety of low-density ionized regions. The predicted 
emission-line ratios can be compared directly with observations of the Galactic WIM (§2), 
while the ionization fractions are useful for analyzing absorption-line data. 

Several caveats regarding these models are worth mentioning. First, the atomic data 
incorporated in CLOUDY, while extensive, is incomplete. For example, the low-temperature 
dielectronic recombination coefficients for most elements in the third and fourth row of 
the periodic table are estimated (Ferland 1996; Ferland et al. 1998). For some elements 
this recombination pathway can be a significant contributor to the ionization balance at 
nebular temperatures. Although the atomic data are somewhat incomplete for many of 
the elements considered, we believe it is important to include the current best values of 
the ionization fractions of all the elements that might be of interest. Second, some of the 
emission-line predictions used to match the WIM observations are highly sensitive to the 
predicted temperature of the medium (Osterbrock 1989). We discuss this further in §3.3. 
CLOUDY calculates the thermal equilibrium structure of a model nebula, but there may be 
additional heating and cooling processes in the WIM that are not taken into account in the 
models. We assume that the WIM is photoionized purely by stellar continuum radiation (as 
well as diffuse radiation from the nebular gas itself). Other processes may play a secondary 
role in ionizing the WIM, such as magnetic reconnection, dissipation of turbulence, or the 



"'^AU emission-line ratios in this paper are intensity ratios, where the intensities are measured in energy 
tinits, not photon units. 
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mixing of cool material with hot, coUisionally ionized gas (Raymond 1992; Slavin, ShuU, & 
Begelman 1993; Minter & Balser 1997). 

3.2. The Standard Model 

The parameter space explored by our models spans a wide range of physical properties 
in the ionized gas. Our preferred model for the WIM, hereafter referred to as the "standard 
model," assumes WNM-hke abundances with interstellar grains included. The ionizing 

spcctrTim is that of a 35,000 K star, consistent with the upper limits to the He I A5876 line 
emission set by observations (Reynolds & Tufte 1995). We assume a very dilute radiation 
field characterized by log(g) = —4.0. 

Table 3 contains the emission-hne intensities relative to Ha predicted by our CLOUDY 
models. We give values for the Xgdge — 0.10, x^dge = 0.95, and standard composite models. 
The composite model provides a reasonable description of the observed emission-line ratios, 
with the possible exceptions of [O I] A6300/Hq; and [O III] A5007/H«. The standard model 
produces less [O III] emission relative to Ha than is suggested by the observations. This is 
a common feature of our models. The predicted [O I]/Hq; ratio may also be slightly low 
compared to the observed values (see Table 1). While most directions have observable 
values of [O Ij/Ho; ~ 0.02, our composite models typically predict ^ 0.015. The strength 
of the [O l] A6300 line, however, is strongly dependent on the temperature of the medium 
(Reynolds et al. 1998a). 

We are not overly concerned about the possible under-production of [O III] A5007 
relative to Ha in our models compared to the observed upper limit of ^ 0.1. This 
observational hmit was derived by noting that the current observed values (Table 1) are 

at best 2-3o" detections, so that the actual value may be considerably lower than this 
limit. Enhanced [O III] emission can be produced by alternate mechanisms (e.g., shocks) as 
well as by photoionization, so we expect that the observed values should indeed provide a 
conservative upper bound for the model predictions. 

Previous investigations seeking to model the WIM emission-line properties (e.g., 
Domgorgen & Mathis 1994) predicted higher [O I] A6300 and [O III] A5007 emission than 
our standard model. To disentangle the differences between these earlier works and our 
models, we have calculated CLOUDY models with the abundances, ionizing spectrum, 
and ionization parameter adopted by Domgorgen & Mathis. We compare their predicted 
emission-line ratios with those from our CLOUDY models in Table 4 and find good 
agreement in most cases. The slight remaining discrepancies, particularly those for [O III], 
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are probably due to differences in the atomic data used in the models (Ferland et al. 1998). 

The elements included in the models compared in Table 4 were limited to H, He, N, O, 
Ne, and S. As a result, the cooling is less efficient than in our standard model, which 
incorporates elements up through Zn, and the resulting nebular temperatures are higher 
in the Domgorgen & Mathis models {Tneb ~ 8000 K versus Tneb ~ 6700 K in our standard 
model). An extensive comparison of CLOUDY to other ionization codes can be found in 
Ferland et al. (1994). 

We present volume-averaged ionization fractions for the standard model in Table 5. 
Values of x{X^) are listed for the first three ionization stages of 20 elements. Ionization 
potentials can be found in Table 2, where we have underlined those ionization stages 
that have at least one observable resonance fine in the 912-3000 A bandpass suitable for 
high-resolution ultraviolet absorption-line work. In the standard composite model, hydrogen 
is almost fully ionized, x{H'^) — 0.81. 

3.3. Sensitivity of Emission-Line Ratios to Model Pcirameters 

The emission-line ratios predicted by our models have varying sensitivities to the input 
model parameters (abundance, ionizing spectrum, ionization parameter). In Figure 1 we 
illustrate the dependences of the emission line ratios on the effective temperature of the 
central ionizing source and the abundance pattern adopted for models with log(g) = —4.0. 
The shaded regions indicate the ranges of observed values as summarized in Table 3. We 
have included [O II] A3727 for completeness, though this line has not yet been observed in 
the Galactic WIM. Of the lines studied, [S II] A6716 and [O I] A6300 are the most sensitive 
to the adopted set of abundances. He I A5876 and [O III] A5007 are affected only slightly 
by the chosen abundances, but vary by at least an order of magnitude over the effective 
temperature range plotted. [S II] A6716 and [N II] A6583 are the least sensitive to the 
ionizing spectrum, varying by less than a factor of 2 over the same range. 

The upper limit to the He I/Ha ratio implies that the ionizing spectrum of the WIM 
produces a ratio of He- to H-ionizing photons less than that produced by an O star with 
Teff ^36,000 K (Reynolds & Tufte 1995). The observed [O I] /Ha constraint is barely 
satisfied over this range of stellar effective temperatures, which indicates that there may be 
some additional nebular heating required to increase the predicted [O I] emission. Note, 
however, that [O I] measurements exist for a small number of WIM sight lines, and there 
appears to be a significant variation in the [O I] /Ha ratio between these directions. The 
assumption that these emission line ratios are representative of the WIM as a whole needs 
to be confirmed observationally. 
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The average nebular electron temperature in our standard model is Tneb — 6700 K, 
which is determined by the input model parameters and the detailed physics of heating and 
cooling in the nebula. We tabulate the effects of nebular temperature on the emission-line 
ratios for the standard composite model in Table 6, where we have forced Tneb to be 
uniform throughout our model nebulae at temperatures of 6000, 8000, and 10000 K. The 
recombination line ratio He I/Ha is relatively unaffected by nebular temperature, as 
expected, but the forbidden emission lines vary greatly in strength as a function of T^eb- 
Changes of ~ 500 — 1000 K in Tneb can increase [O I]/Hq; by a factor of ~ 2. In general, 
values of Tneb as high as 8000 K are inconsistent with the observations in Table 3. But, the 
observed increase in [N II] /Ha and [S II] /Ha with height above the Galactic plane suggests 
temperatures as high as T^eb ~ 10^ K may be appropriate at large z-distances (Haffner 
et al. 1999). Observations of [O II] A3727 could provide an excellent test for additional 
heating sources since the [O II] /Ha ratio is expected to vary by a factor of about 16 for 
6000 ^T„e6^ 10000 K. 

We note that grain heating in our models is approximately balanced by grain cooling 
under the conditions considered. The heating rate due to ejection of photoelectrons from 
grains is less than 10% of the H I bound-free photoionization heating in our models. 
We refer the reader to Cox & Reynolds (1992) and Baldwin et al. (1991) for additional 
discussions of grain heating in the ionized ISM. 

In general, the emission-line ratios do not depend strongly upon the value of q adopted 
for the standard model. A value of log(g) = —3.5 versus —4.0 results in changes of less than 
10% for [N II]/Ha, [S II]/Ha, [O I]/Ha, and He I/Ha. The [O III]/Ha ratio varies by less 
than ~40%, with smaller changes occurring for values of Tgff < 38, 000 K. 

3.4. Sensitivity of Ionization Fractions to Model Pcirameters 

The ionization fractions for the standard model (Table 5) have modest dependencies 
on the input parameters used to model the emission line observations. We illustrate these 
dependencies in Figures 2 and 3, where we plot the ionization fractions as functions of Te// 
and log(g), respectively. The standard Xedge = 0.10, Xedge = 0.95, and composite models 
are shown with dashed, dotted, and solid lines, respectively. The ionization stages are 
differentiated by the type of symbol used, with open circles for neutral species, filled circles 
for singly ionized species, and crossed circles for doubly ionized species. For a few ions, the 
ionization fractions depend strongly upon Tg// (e.g., Hel-II, C III), but in most other cases 
of interest here, the variations with T^ff are modest (Aloga;(X*) < 0.2 dex). Similarly, 
changes in the adopted reference abundances do not change the ionization fractions 
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significantly. The variations of logx{X'^) with log{q) shown in Figure 3 are typically less 
than 0.2 dex for values of log(g) between -4.0 and -3.0. 

The ionization fractions do not depend significantly on the nebular temperature since 
their values are governed in large part by recombination processes. Changes in x{X^) are 
less than ~ 10% for Tneb — 6,000 — 10, 000 K. The lack of variation in the strength of 
the permitted He I A5876 recombination line relative to Ha with nebular temperature is 
another manifestation of this insensitivity (see Table 6). 

The information provided in Table 5 and Figures 2 and 3 should have a broad range 
of applicabihty to studies of the WNM and WIM. The ionization fractions can be scaled 
easily to other preferred values of the input model parameters. As better emission-line 
data become available and additional constraints are set on the properties of the WIM, 
the standard model ionization fractions can also be refined using the information presented 



4. Application to Absorption-Line Abundance Studies of the WNM and WIM 



In this section we outline a simple strategy for estimating elemental abundances in 
the warm neutral and ionized media of the Galaxy. For simplicity, any cold neutral gas 
along the sight line can be considered as WNM material in the following discussion. Many 
singly ionized elements can be found in both the WNM and WIM (e.g., Al II, Si II, P II, 
S II, Fe II), making it necessary to account for the relative contributions of each type of gas 
to the observed column densities. Higher ionization stages are preferentially found in the 
WIM, and these can be used to isolate the elemental abundances in the WIM. 

Assume that an ion is found only in the ionized gas and presumably has no 
appreciable abundance in the WNM (e.g., S+^). Then, we have 



where iV(X*) is the observed column density of X*. In general, the total amount of an ion 
is the sum of contributions from the WIM and WNM. 



here. 



4.1. Methodology 



N{X')wiM = N{X') 

N{X')wNM = 




(3) 



N{X') = N{X')wiM + N{X') 

WNM 



(4) 
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The column density of element X in the WIM is given by 



N{X)wiM = N{X^)/x{X^) (5) 

where x{X^) is the volume-averaged ionization fraction found in Table 5. The ionization 
fractions in Table 5 are appropriate for the average properties of the WIM. The 
corresponding logarithmic abundance of the element relative to hydrogen in the WIM is 



A{X)y^IM = log ^^^^ 



N{H) 



+ 12.00 (6) 



WIM 



It is difficult to estimate N{H)-wim directly since N{H^) is not directly observed, so a 
suitable proxy for N{H)wim may need to be found. 

In the WNM, the column density of element X is equal to the observed value less the 
amount in the WIM, which is simply the sum of the individual contributions of X' to N{X) 
in the WNM. 

N{X)^NM = E (n{X^) - ^fi^n^ (7) 

In Eq.(7), the second term in parentheses accounts for the amount of X* in the WIM. We 
use the following equations to estimate the amount of H I in the WIM (Eq. 8) and WNM 
(Eq. 9): 

N{H'>)wNM = iV(ifO) [l - x{H'')Fh] (9) 

In Eqs. (8) and (9), Fh is the fraction of hydrogen along the sight line associated with the 
WIM. The ionization fraction of H I in the WIM, x{H^), in our models is ~ 5 — 30%, so 
most of the H I is associated with the WNM. 



Combining Eqs. (6) and (8), we derive the logarithmic abundance of X relative to H in 
the WNM: 



+ 12.00 

WNM 



(10) 



4.2. An Example: The HD 93521 Sight Line 

To demonstrate our approach we apply the method outhned in Eqs. (3)-(10) to 
estimate the gas-phase abundance of S in the warm halo clouds towards HD 93521. Spitzer 
& Fitzpatrick (1993) have presented high-resolution observations of this sight line obtained 
with the Goddard High Resolution Spectrograph on-board the HST. The neutral and 
ionized species occur at very similar velocities along this sight line, so the contributions of 
ionized gas to the derived elemental abundances may be important. Spitzer & Fitzpatrick 
(1993) measured the column densities of S II, S III, and H I in the individual "warm fast" 
clouds along this sight line. In the first five columns of Table 7 we give their identifications, 
LSR velocities, and column densities of S II, S III, and H I for each of these clouds. These 
values are followed in column 6 by the column density of S II in the WIM calculated using 
Eq. (7) and the model predictions given in Table 5. The ratio of S II in the WNM to total 
S II in each cloud is listed in column 7. Equations (8) and (9) were used to derive the 
ratio of H I in the WNM to the observed value. We list this ratio in column 8. Column 
9 contains the derived abundance of S in the WNM found from Eq. (10), and column 10 
contains A{S) — Ac{S), the logarithmic gas-phase abundance of sulfur relative to the WNM 
reference abundance given in Table 2. 

Table 7 shows the results of these calculations using both the standard composite 
model and the Xedge — 0.10 model results for x{X^). The cloud-to-cloud scatter in the values 
of A{S) — Ac{S) is generally of the order of the measurement errors on A'"(S III). Values of 
A{S) — Ac{S) for components 1,2, and 4 for the standard composite model are consistent 
with a negligible gas-phase depletion of S onto dust grains. The value for component 3 
suggests either that there may be a small amount of S incorporated into dust or that 
the ionization properties or intrinsic abundances of cloud 3 vary enough to produce the 
observed deviation from the reference abundance. In the Xedge = 0.10 model, the -0.14 dex 
deviation from the reference value is a 2a excursion. The cloud 3 value is 0.2 dex below the 
value of A{S) derived by Spitzer & Fitzpatrick (1993), who assumed that all of the S II in 
these clouds is associated with the H I (i.e., they did not apply any ionization corrections). 

We find that a substantial ionization correction is necessary for S II to explain the 
amount of S III observed along the sight line. We find that roughly 50% of the S II is in 
the WNM, with the other 50% contained in the WIM (see Table 7). It is interesting that 
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Spitzer & Fitzpatrick (1993) found an electron density for the HD 93521 clouds comparable 
to that of the WIM. The main difference in addressing ionization in these two studies is 
that they used the observed ratio N(S II)/N(S III) and measured electron density to derive 
the photoionization rate for S II under conditions of photoionization equilibrium, whereas 
our results take into account the explicit ionizing spectrum required to match the WIM 
properties. We also predict varying ionization conditions between clouds, whereas they 
assumed constant properties. Spitzer & Fitzpatrick (1993) argued that the free electrons 
and neutral gas along the HD 93521 sight hne are well mixed. Our calculation is consistent 
with this result, but does not require it. Additional observations of other ionized gas 
species in the HD 93521 clouds would help to characterize their ionization conditions and 
the impact of ionization on the derived elemental abundances. 

Analyses such as these have also been used successfully to study the abundances in 
the WNM of the Galactic disk and halo and in intermediate- velocity clouds (e.g., Sembach 
1995; Howk et al. 1999). Eventually, similar studies of the WIM abundances could be made 
(see Howk & Savage 1999). Wc note that especially good indicators of partially ionized 
gas are observable at the far- ultraviolet wavelengths to be observed by FUSE. Information 
about spectral lines that are useful probes of photoionized gas can be found in Sembach 
(1999) and Howk & Savage (1999) - see Table 1 in both papers. 



5. Concluding Remeirks 

The outlook for a better understanding of the abundances and physical properties of 
the WIM is bright. Advances in measuring the emission-line properties of the WIM will 
probably occur at optical and near-infrared wavelengths in the near future. Emission-line 
studies of the WIM at wavelengths A < 3000 A will require sensitive detectors with large 
fields-of-view. The predicted intensities of the strongest emission lines at far-ultraviolet 
wavelengths (A < 1200 A) in our models are currently below the sensitivity of FUSE. In the 
standard composite model [C II] A1020, C III A977, N II A1085, and N III A990 all have 
intensities less than 10~'^ Ina- At these strengths, the far- ultraviolet WIM emission lines 
measured by FUSE in its 30"x30" aperture should fall well below the expected background 
levels of ~ 0.5 counts scc^^ cm^^ and should pose no significant contamination problem for 
studies of hot gas emission (e.g., O VI) from the highly ionized ISM. Absorption line studies 
of the WIM with FUSE should be feasible, with excellent diagnostics of photoionized gas 
ranging in ionization energies from 1 to 4 Rydbergs (e.g., Cll-lll, Nl-lll, Arl-ll, Fell-lll, 
Pll-lV, Slll-iv). 

The results of this work are applicable to warm ionized gas regions in other galaxies. A 
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companion study of warm ionized gas and its effects on elemental abundance determinations 
in damped Lyman-a absorbers can be found in Howk & Sembach (1999). 

We thank Gary Ferland and his co-workers at the University of Kentucky for all of 
their work on the CLOUDY code. We also thank John Mathis for useful discussions about 
modeling the WIM and Ron Reynolds for helpful suggestions for improving the manuscript. 
KRS and JCH acknowledge support from NASA Long Term Space Astrophysics grant 
NAG5-3485 and grant GO-07270.01-96A from the Space Telescope Science Institute, which 
is operated by the Association of Universities for Research in Astronomy, Inc., under NASA 
contract NAS5-26555. 



A. Ionization Corrections for Low-Density H ii Regions 

Although the main purpose of this work is to determine ionization fractions for the 
WIM of the Galaxy, our CLOUDY calculations also have applicability to studies of ionized 
gas in low-density H II regions. H II regions can be a significant source of singly-ionized 
species, which are also the primary tracers of many elements in the neutral ISM. Howk et 
al. (1999) have discussed the effects of H II region gas on the the study of the gas-phase 
abundances in the WNM along the // Columbae sight line. 

The data given in Table 5 describe the volume-averaged properties of our model 
nebulae. For a sight line through a diffuse H II region toward the central ionizing 
source, it may be more appropriate to use radially-averaged ionization fractions (see 
Howk & Savage 1999). In Table 8 we hst the radially- averaged ionization fractions for 
all 20 elements considered in this study. We provide ionization fractions for diffuse H II 
regions characterized by log (g) = —3.0 and —4.0, with an ionizing spectrum from an 
ATLAS model atmosphere with T^ff = 35,000 K. We present model results for both 
Xedge = 0.10 and 0.95. Figures 4a and 4b show the variations in these ionization fractions 
as a function of input stellar effective temperature for a select number of elements for the 
log (q) = —4.0 models. For the elements Al, Si, Mg, Cr, and Mn, the variation in the 
ionization fractions of the singly and doubly ionized stages is less than 0.05 dex over the 
range 31, 000 < Te// < 41, 000 K. Figure 5 shows the variation in ionization fraction with 
log (g), assuming an input spectrum with T^ff — 35, 000 K. 

Howk & Savage (1999) have tabulated radially- average ionization fractions for a limited 
number of ions and discussed in detail how the ionization fractions vary with assumed 
model parameters. The results presented below are more thorough in the sense that more 
elements are covered, though over a smaller range of ionization parameter. These ionization 
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fractions are appropriate for very low-density H II regions. In particular, these results 
should be useful for studying contamination from H II regions surroTinding early-type stars 
at large distances from the Galactic plane or along short, low-density sight lines where an 
appreciable number of additional H II regions has not been intercepted. 

The application of the ionization fractions given in Table 8 should be approached 
in a method similar to that described in §4. When possible, multiple observations of 
adjacent ions should be used to constrain the most appropriate model parameters. The 
radially-averaged values tend to favor highcr-ionization gas than the volume-averaged 
results, since the partially-ionized outer regions of the nebula are given less weight in the 
averaging process. 
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Figure Captions 

Fig. 1. — Predicted emission line ratios versus the assumed effective temperature of tlie 
underlying stellar ionizing source from the CLOUDY model nebula calculations. The three 
curves in each panel depict the relationships for the different abundance patterns indicated 
in the legend in the lower left panel. For each set of reference abundances we show the 
composite model with an ionization parameter log{q) — —4.0. Shaded regions indicate the 
range of observed emission-line ratios (see Table 1). Note that [O II] A3727 emission has not 
yet been observed in the WIM. 

Fig. 2. — Volume-averaged ionization fractions, logx{X'), versus input stellar effective 
temperature of the ionizing central source in our nebular models with WNM abundance 
patterns. The curves shown are for the Xedge — 0.10 (dashed hnes), Xedge — 0.95 (dotted 
hnes), and composite (sohd lines) cases. The ionization stage (I-III) is indicated by the type 
of symbol plotted: open circles (I), filled circles (II), and crossed circles (III), 
a) He, C, N, and O; b) S, P, Ar, and Zn; c) Si, Mg, Mn, and Cr; d) Fe, Ti, Al, and Ni. 

Fig. 3. — Volume-averaged ionization fractions, loga;(X''), versus ionization parameter, 
\og{q), for a model nebula with a WNM abundance pattern and an ionizing source 
characterized by Te// = 35,000 K. The curves shown are for the Xedge — 0.10 (dashed 
lines), Xedge = 0.95 (dotted lines), and composite (solid lines) cases. The ionization stage 
(I-III) is indicated by the type of symbol plotted: open circles (I), filled circles (II), and 
crossed circles (III). 

Fig. 4. — Radially- averaged ionization fractions, logx{X^), versus effective temperature 
assumed for the ionizing central star. The curves shown are for the Xedge — 0.10 (dashed 
lines), Xedge = 0.95 (dotted lines) models. All models assume WNM abundance patterns, but 
the ionization fractions are relatively insensitive to the adopted abundances. The ionization 
stage (I-III) is indicated by the type of symbol plotted: open circles (I), filled circles (II), 
and crossed circles (III), 
a) C, N, O, and Zn; b) S, P, Fe, and Al. 

Fig. 5. — Radially-averaged ionization fractions, loga:(X*), versus ionization parameter, 
log(g). The curves shown are for the Xedge = 0.10 (dashed lines), Xedge = 0.95 (dotted lines) 
models. These models assume an ionizing central star with Teff = 35,000 K. All models 
assume WNM abundance patterns, but the ionization fractions arc relatively insensitive to 
the adopted abundances. The ionization stage (I-III) is indicated by the type of symbol 
plotted: open circles (I), filled circles (II), and crossed circles (III). 
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TABLE 1 

Emmission Line Measurements of the Warm Ionized Medium 



Direction 
(/, b) 


[Sii] A6716 
Ha 


[ATii] A6583 
Ha 


[Oi] A0300 
Ha 


[Oiii] A5007 
Ha 


Hei A5876 
Ha 


Reference'^ 


96?0, 0?0 


0.29±0.05 


0.28±0.06 




0.058±0.025 




1,2 


114?0, 0?0 


0.24±0.04 




0.020±0.003 




<0.02 


2,3,4 


130?0, 0?0 


0.27±0.06 




0.028±0.009 






2,3 


184?4,-0?6 1 
186?0,0?0 


0.23±0.04 


0.21±0.02 




<0.06 




2,5 


194?0, 0?0 


0.38±0.05 


0.44±0.05 




0.068±0.025 


<0.02 


1,2,4 


97?3,-29?7 1 
97?7,-26?7 


0.39±0.11 






<0.20 




1,2 


144?0, -21?0 


0.44 


0.70 


< 0.12 


<0.57 




6 



(weak shock) 

123 to 164, -6 to -35 0.2-0.4'= 0.4-0.6'' 

(Perseus Arm) 



< 0.012^= 
0.044 ±0.011' 



3,7 



^References: (1) Reynolds 1985b, (2) Reynolds 1985a, (3) Reynolds et al. 1998a, (4) Reynolds & Tufte 1995, (5) 
Reynolds, Roessler, & Scherb 1977; (6) Ogden & Reynolds (1985), (7) Haffner, Reynolds, & Tufte 1999. 

''The ranges of values quoted are appropriate for Perseus Arm gas within <^ 1 kpc of the Galactic plane having an Ha 
intensity of ~ 2 — 3 R. There is a slightly larger scatter to higher values in the Local Arm gas in the same directions. 
Both the [S ii]/Ha and [N ii]/Ha ratios are anti-correlated with Iho- The upper envelope of observed ratios extends to 
values several times higher than those listed. 

■^These values are for two velocities (-60 and -31 km s^^) at the position / — 130°, b — —7.5°. 
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TABLE 2 

Ionization Potentials and Reference Abundances 



Element Ionization Potentials (eV)** A{X)^ 





T 
1 


TT 


TTT 
111 


ooiai 




VV IN iVl 


H 


13.60 






12.00 


12.00 


12.00 


He 


24.59 


54.42 




11.00 


11.00 


11.00 


C 


11.26 


24.38 


47.89 


8.55 


8.27 


8.15 


N 


14.53 


29.60 


47.45 


7.97 


7.80 


7.88 





13.62 


35.12 


54.93 


8.87 


8.56 


8.50 


Ne 


21.56 


40.96 


63.45 


8.07 


8.20 


8.07 


Na 


5.14 


47.29 


71.64 


6.31 


6.11+ 


6.11 


Mg 


7.65 


15.04 


80.14 


7.58 


7.38 


7.16 


Al 


5.99 


18.83 


28.45 


6.48 


6.19 


5.18 


Si 


8.15 


16.34 


33.49 


7.55 


7.20 


7.18 


P 


10.49 


19.72 


30.18 


5.57 


5.37'f 


5.37 


s 


10.36 


23.33 


34.83 


7.27 


6.97 


7.07 


Ar 


15.76 


27.63 


40.74 


6.52 


6.32+ 


6.32 


Ca 


6.11 


11.87 


50.91 


6.36 


6.16+ 


5.06 


Ti 


6.82 


13.58 


27.49 


4.93 


4.73+ 


3.86 


Cr 


6.77 


16.50 


30.96 


5.68 


5.48+ 


4.55 


Mn 


7.43 


15.64 


33.67 


5.53 


5.33+ 


4.56 


Fe 


7.87 


16.18 


30.65 


7.51 


7.40 


6.21 


Ni 


7.64 


18.17 


35.17 


6.25 


6.05+ 


4.95 


Zn 


9.39 


17.96 


39.72 


4.65 


4.45+ 


4.45 



Ionization potentials are from Moore (1970). Underlined values indicate that the ion- 
ization stage has at least one ground-state resonance line in the FUSE/HST ultraviolet 
bandpass (912-3000 A) that could be observable in the WIM. 

''Reference abundances are logarithmic values on a scale where A{X) = \og{X/H) + 12.00. 
Solar abundances arc metcoritic values from Anders & Grevesse (1989), except for C, N, 
and O, which are photospheric values from Grevesse & Noels (1993). B-star abundances 
are from Kilian-Montenbruck et al. (1994), except for those cases marked with a f, which 
are equal to the solar value minus 0.2 dex. The WNM (warm neutral medium) abundances 
are equal to the solar abundances modified by dust as observed in the warm diffuse clouds 
toward jj, Columbac (Howk ct al. 1999). WNM values for C, N, and O are from Sofia et al. 
(1997), Meyer et al. (1997), and Meyer et al. (1998), respectively. 
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TABLE 3 

Standard Model WIM Emission Lines 



Line 


Volume- Averaged Intensity Ratio (7^) 


Observed** 


Xedge = 0.10 


Xedge = 0.95 


Composite'^ 


[S ii] A6716 


0.258 


0.325 


0.292 


0.2-0.4 


[N ii] A6583 


0.439 


0.462 


0.451 


0.3-0.6 


[0 i] A6300 


0.005 


0.021 


0.013 


0.01-0.03 


[0 III] A5007 


0.002 


0.002 


0.002 


^0.1 


He I A5876 


0.017 


0.013 


0.015 


<0.02 



'^This summary of observed intensity ratios is based on the data listed 
in Table 1. 

^The composite model is comprised of equal amounts of Xgdge = 0.10 
and Xedge = 0.95 material. 
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TABLE 4 
Emission Line Comparison 



Line 


Intensity Ratio' 




CLOUDY 


DM94 


[S ii] A6716 


0.77 


0.66 


[N it] A6583 


0.37 


0.37 


[0 i] A6300 


0.051 


0.055 


[0 ii] A3727 


1.17 


1.18 


[0 III] A5007 


0.009 


0.038 


He I A5876 


0.025 


0.029 




8040 


8200 



'^These volume-averaged intensity ratios 
are predictions from the CLOUDY and 
Domgorgcn & Mathis (1994) ionization codes 
for a model nebula with Xedge = 0.95, T^e// 
= 38,000 K, and log(g) = -4.0. The calcula- 
tions included H, He, N, O, Ne, and S scaled 
to the Orion abundances from Peimbert et al. 
(1992). 

'^Volume averaged nebular temperature de- 
fined by Tneb = J TuedV/ J UedV. 



24 



TABLE 5 

Ionization Fractions for the Standard Model**'^ 



Element Xedge = 0.10 Xedge = 0.95 Composite'^ 





I 


II 


III 


I 


II 


III 


I 


II 


III 


H 


-1.30 


0.02 




0.49 


0.17 




0.73 


0.09 




He 


-0.17 


-0.50 




-0.07 


-0.83 




-0.12 


-0.63 




C 


-2.24 


-0.02 


-1.31 


-2.16 


-0.01 


-1.63 


-2.20 


-0.02 


-1.44 


N 


-1.33 


-0.03 


-1.72 


-0.52 


-0.16 


-2.05 


-0.76 


-0.09 


-1.85 





-1.26 


-0.03 


-2.74 


-0.48 


-0.18 


-3.07 


-0.71 


-0.09 


-2.87 


Ne 


-0.97 


-0.05 


-3.88 


-0.54 


-0.15 


-4.21 


-0.70 


-0.10 


-4.02 


Na 


-1.93 


-0.01 


-5.53 


1.75 


0.01 


5.82 


1.83 


0.01 


-5.65 


Mg 


-2.30 


-0.15 


-0.54 


-2.22 


-0.07 


-0.82 


-2.26 


-0.11 


-0.66 


Al 


-2.81 


-0.03 


-1.16 


-2.83 


-0.02 


-1.44 


-2.82 


-0.03 


-1.28 


Si 


-3.63 


-0.06 


-0.92 


-3.51 


-0.03 


-1.23 


-3.57 


-0.04 


-1.05 


P 


-2.39 


-0.09 


-0.72 


-2.36 


-0.05 


-1.01 


-2.37 


-0.07 


-0.84 


S 


-2.83 


-0.11 


-0.65 


-2.76 


-0.06 


-0.90 


-2.80 


-0.09 


-0.75 


Ar 


-1.86 


0.12 


0.63 


0.94 


0.11 


-0.96 


1.19 


-0.12 


-0.77 


Ca 


-2.33 


-0.71 


-0.10 


-2.31 


-0.65 


-0.11 


-2.32 


-0.68 


-0.10 


Ti 


-3.19 


-0.34 


-0.28 


-2.87 


-0.15 


-0.54 


-3.00 


-0.23 


-0.39 


Cr 


-3.42 


-0.41 


-0.21 


-3.12 


-0.17 


-0.48 


-3.24 


-0.28 


-0.33 


Mn 


-1.68 


-0.31 


-0.31 


-1.39 


-0.16 


-0.58 


-1.51 


-0.23 


-0.43 


Fe 


-3.41 


-0.47 


-0.19 


-3.02 


-0.19 


-0.45 


-3.17 


-0.31 


-0.30 


Ni 


-1.98 


-0.11 


-0.68 


-1.89 


-0.05 


-1.00 


-1.93 


-0.08 


-0.81 


Zn 


-2.15 


-0.04 


-1.08 


-1.96 


-0.02 


-1.40 


-2.05 


-0.03 


-1.21 



'^Thc ionization fraction of an ion is given by x{X^) = N (X^) / N (XTotai) ■ We 
present volume-averaged values of logx(X*) for the first three ionization states of each 
element. 

^The standard model assumes WNM abundances (see Table 1), \og(q) = -4.0, and a 
stellar input spectrum characterized by Tgjf = 35,000 K and log (g) = 4.0. 
■^The composite model is described by 50%{xedge = 0.10) + 50%{xedge = 0.95). 
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TABLE 6 

Effects of Nebular Temperature on Emission Lines 



Line 




Intensity Ratio^ 




(i.OOOK 


S.OOOK 


10. ()()() K 


[S ii] A6716 


0.200 


0.539 


0.990 


[N ii] A6583 


0.291 


0.820 


1.571 


[0 i] A6300 


0.009 


0.032 


0.075 


[0 III] A5007 


0.001 


0.002 


0.005 


He I A5876 


0.016 


0.015 


0.014 


[0 ii] A3727 


0.186 


1.050 


3.031 


[S III] A9532 


0.085 


0.173 


0.271 



^These volume-averaged intensity ratios are for the standard 
composite model, which is comprised of equal amounts of Xedge = 
0.10 and Xedge = 0.95 material. 
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TABLE 7 

Standard Model Results for the HD 93521 Sight Line'' 



Comp. <Vlsr> log7V(Sii) logiV(Siii) logiV(Hi) log7V(Sii)^i„ ^^^wwf^ A(S)-Ac(S) 

(km s-i) 



Composite Model 



1 -66.3 14.09±0.03 13.05±0.21 18.85 13.71 0.58 0.92 7.04 -0.03 

2 -57.8 14.50±0.02 13.11±0.22 19.26 13.77 0.81 0.97 7.17 +0.10 

3 -51.2 14.53±0.03 13.71±0.07 19.29 14.37 0.31 0.87 6.79 -0.28 

4 -38.8 14.02±0.07 13.14±0.31 18.78 13.80 0.40 0.89 6.89 -0.18 
Sum ••• 14.95±0.02 13.95±0.08 19.71 14.61 0.54 0.91 7.01 -0.06 



Xedge = 0.10 

1 -66.3 14.09±0.03 13.05±0.21 18.85 13.59 0.68 0.94 7.08 +0.01 

2 -57.8 14.50+0.02 13.11+0.22 19.26 13.65 0.86 0.97 7.18 +0.11 

3 -51.2 14.53+0.03 13.71+0.07 19.29 14.25 0.48 0.90 6.93 -0.14 

4 -38.8 14.02+0.07 13.14+0.31 18.78 13.68 0.54 0.91 6.98 -0.09 
Sum ••• 14.95+0.02 13.95+0.08 19.71 14.49 0.65 0.94 7.06 -0.01 



^Observed quantities (< Vlsr >, N{S ii), A^(S iii), iV(H :)) are measurements for the "warm, fast" halo clouds observed 
toward HD 93521 by Spitzer & Fitzpatrick (1993). The remaining derived quantities are defined in the text. 
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TABLE 8 

Ionization Fractions for Low-Density H ii Regions'* 



Element 






log (g) = 


-3.0 










log [q) = 


-4.0 






Xedge = 0.10 


X, 


edge = 0.95 




Ige = 0.10 




Ige = 0.95 


I 


II 


III 


I 


II 


III 


I 


II 


III 


I 


II 


III 


H 


-1.81 


-0.01 




-1.08 


-0.04 




-1.60 


-0.01 




-0.87 


-0.06 




He 


-0.44 


-0.20 


-8.06 


-0.34 


-0.26 


-8.13 


-0.42 


-0.21 


-7.52 


-0.28 


-0.32 


-7.63 


C 


-2.81 


-0.17 


-0.50 


-2.71 


-0.14 


-0.57 


-2.52 


-0.12 


-0.61 


-2.37 


-0.10 


-0.72 


N 


-1.87 


-0.14 


-0.59 


-1.13 


-0.15 


-0.66 


-1.65 


-0.10 


-0.72 


-0.90 


-0.14 


-0.83 





-1.78 


-0.05 


-1.04 


-1.07 


-0.08 


-1.11 


-1.57 


-0.04 


-1.18 


-0.85 


-0.09 


-1.30 


Ne 


-1.49 


-0.02 


-1.79 


-1.11 


-0.04 


-1.86 


-1.26 


-0.03 


-1.91 


-0.86 


-0.07 


-2.02 


Na 


-2.28 


-0.00 


-3.02 


-2.26 


-0.00 


-3.08 


-2.14 


-0.00 


-2.88 


-2.01 


-0.01 


-2.99 


Mg 


-2.71 


-0.45 


-0.19 


-2.63 


-0.36 


-0.25 


-2.52 


-0.35 


-0.26 


-2.44 


-0.25 


-0.36 


Al 


-3.27 


-0.18 


-0.58 


-3.29 


-0.15 


-0.65 


-3.05 


-0.13 


-0.70 


-3.03 


-0.10 


-0.80 


Si 


-4.05 


-0.30 


-0.32 


-4.04 


-0.24 


-0.39 


-3.88 


-0.21 


-0.42 


-3.78 


-0.16 


-0.53 


P 


-3.02 


-0.34 


-0.29 


-2.94 


-0.28 


-0.35 


-2.67 


-0.26 


-0.37 


-2.56 


-0.19 


-0.48 


s 


-3.48 


-0.38 


-0.24 


-3.37 


-0.31 


-0.30 


-3.13 


-0.28 


-0.33 


-2.98 


-0.21 


-0.43 


Ar 


-2.42 


-0.36 


-0.26 


-1.57 


-0.30 


-0.32 


-2.21 


-0.33 


-0.28 


-1.33 


-0.26 


-0.40 


Ca 


-3.08 


-1.21 


-0.03 


-3.05 


-1.14 


-0.03 


-2.64 


-0.97 


-0.05 


-2.56 


-0.85 


-0.07 


Ti 


-3.74 


-0.74 


-0.20 


-3.54 


-0.55 


-0.25 


-3.51 


-0.59 


-0.22 


-3.21 


-0.39 


-0.31 


Cr 


-4.04 


-0.89 


-0.12 


-3.81 


-0.63 


0.17 


3.76 


-0.70 


0.14 


3.46 


0.44 


0.23 


Mn 


-2.32 


-0.67 


-0.15 


-2.05 


-0.51 


0.20 


-2.00 


-0.53 


-0.19 


-1.69 


-0.37 


-0.28 


Fe 


-4.03 


-0.90 


-0.16 


-3.77 


-0.64 


-0.21 


-3.75 


-0.73 


-0.17 


-3.38 


-0.47 


-0.26 


Ni 


-2.56 


-0.47 


0.21 


-2.51 


-0.37 


0.27 


2.28 


0.34 


0.30 


-2.17 


-0.24 


-0.41 


Zn 


-2.57 


-0.22 


-0.41 


-2.54 


-0.18 


0.48 


2.41 


0.17 


0.50 


-2.26 


-0.13 


-0.61 



''The ionization fraction of an ion X* is given by x{X^) = N{X'^)/N{XTotsd)- We present radially-averaged 
values of logx{X'^) for the first three ionization states of each element. Each model assumes WNM abundances (see 
Table 2) and a stellar input spectrum characterized by Te// = 35,000 K and log(g) = 4.0. Results are listed for 
ionization parameters of log (g) = -3.0 and -4.0. 
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